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Abstract 

A  newly  developed  microstructure-sensitive  extreme  value  probabilistic  framework  for  fatigue  variability  based 
on  computational  polycrystal  plasticity  is  exercised  to  compare  the  driving  forces  for  fatigue  crack  formation 
(nucleation  and  early  growth)  at  room  temperature  for  four  different  microstructure  variants  of  duplex  a  Ti-6AI-4V 
alloy.  The  aforementioned  probabilistic  framework  links  certain  extreme  value  fatigue  response  parameters  with 
microstructure  attributes  at  fatigue  critical  sites  through  use  of  marked  correlation  functions.  By  applying  this 
framework  to  study  the  driving  forces  for  fatigue  crack  formation  in  these  microstructure  variants  of  Ti-6AI-4V, 
these  microstructures  can  be  ranked  in  terms  of  relative  high  cycle  fatigue  (HCF)  performance  and  the  correlated 
microstructure  attributes  that  have  the  most  influence  on  the  predicted  fatigue  response  can  be  identified. 
Nonlocal  fatigue  indicator  parameters  (FIPs)  based  on  the  cyclic  plastic  strain  averaged  over  domains  on  the  length 
scale  of  the  microstructure  attributes  (e.g.  grains,  phases)  are  used  to  estimate  the  driving  force(s)  for  fatigue  crack 
formation  at  the  grain  scale.  By  simulating  multiple  statistical  volume  elements  (SVEs)  using  crystal  plasticity 
constitutive  relations,  extreme  value  distributions  of  the  predicted  driving  forces  for  fatigue  crack  formation  are 
estimated  using  these  FIPs.  This  strategy  of  using  multiple  SVEs  contrasts  with  simulation  based  on  a  single 
representative  volume  element  (RVE),  which  is  often  untenably  large  when  considering  extreme  value  responses. 
The  simulations  demonstrate  that  microstructures  with  smaller  relative  primary  a  grain  sizes  and  lower  volume 
fractions  of  the  primary  a  grains  tend  to  exhibit  less  variability  and  smaller  magnitudes  of  the  driving  forces  for 
fatigue  crack  formation.  Additionally,  the  extreme  value  FIPs  are  predicted  to  most  likely  occur  at  clusters  of 
primary  a  grains  oriented  for  easy  basal  slip.  In  some  cases,  surrounding  grains/phases  with  soft  orientation  shed 
load  to  less  favorably  oriented  primary  a  grains,  producing  extreme  value  FIPs. 

Keywords:  high  cycle  fatigue,  extreme  value  statistics,  Ti-6AI-4V,  microstructure,  crystal  plasticity 
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1.  Introduction 

Variability  of  the  fatigue  life  of  engineered  components  arises  from  microstructure  stochasticity.  This  is 
particularly  true  for  the  processes  of  fatigue  crack  nucleation  and  microstructurally  small  crack  growth  in  ductile 
metallic  material  systems,  which  are  driven  by  localized  plasticity  and  the  accumulation  of  dislocations  against 
obstacles  or  the  development  of  particular  dislocation  structures  (e.g.,  persistent  slip  bands)  [1].  In  high  cycle 
fatigue  (HCF)  of  metals  for  which  cyclic  stress  amplitudes  are  below  the  macroscopic  yield  stress,  plasticity  is  quite 
heterogeneous  and  localized  at  microstructure  attributes  that  raise  the  local  stresses  sufficiently  to  induce  flow  [2], 
are  favorably  oriented  for  preferential  yield,  or  both.  Specific  life  limiting  attributes  can  vary  from  specimen  to 
specimen  or  from  component  to  component  due  to  the  stochasticity  of  material  microstructure.  Thus,  scatter  in 
the  HCF  life  of  specimens  or  components  depends  on  the  extreme  value  probabilities  of  existing  microstructure 
attributes  that  elevate  the  driving  forces  for  fatigue  crack  formation,  constituting  size  effects. 

1.1  Extreme  value  statistics  and  metal  fatigue 

Current  methods  that  estimate  fatigue  resistance  as  a  function  of  the  extreme  value  statistics  of  microstructure 
are  limited  to  distributions  of  a  single  microstructure  attribute;  these  methods  do  not  consider  the  effect  of 
multiple  interacting  microstructure  attributes  on  the  extreme  value  response  distributions  (e.g.,  fatigue  crack 
formation).  For  example,  Atkinson  and  Shi  [3]  reviewed  various  extreme  value  models  that  predict  the  fatigue 
resistance  in  clean  steels  based  on  the  largest  inclusions,  which  control  fatigue  crack  formation  in  these  materials. 
The  first  class  of  models  discussed  is  based  on  the  log-normal  distribution  (or  any  standard  distribution  that 
considers  the  entire  population)  of  inclusion  size.  Accuracy  of  these  models,  however,  is  limited  by  the 
experimental  difficulty  of  correctly  capturing  the  tails  of  the  distribution  for  inclusion  size  because  of  the  paucity  of 
extreme  value  data.  The  second  class  of  models  discussed  by  Atkinson  and  Shi  is  based  on  the  classical  extreme 
value  Gumbel  [4]  distribution.  For  example,  Beretta  and  Murakami  [5]  used  classical  Gumbel  extreme  value 
statistics  to  estimate  the  size  of  the  largest  inclusion  based  on  specimen  volume,  which  was  subsequently  used  to 
estimate  fatigue  strength.  The  third  class  of  models  is  based  on  the  generalized  Pareto  distribution.  Atkinson  and 
Shi  argue  that  the  generalized  Pareto  distribution  is  more  effective  than  the  two  previously  described  methods 
because  it  incorporates  a  limit  on  the  maximum  inclusion  size,  whereas  the  extreme  value  Gumbel  distribution 
predicts  a  monotonically  increasing  inclusion  size  for  increasing  volume  of  material  sampled.  However,  it  is  very 
difficult  to  estimate  the  upper  bound  on  this  type  of  attribute.  Additionally,  all  three  model  classes  do  not  consider 
how  interactions  between  the  inclusion  and  surrounding  matrix  or  microstructure  affect  fatigue  performance. 
Although  first  order  approaches  based  on  a  single  microstructure  attribute  (e.g.,  inclusion  size)  may  be  sufficient 
for  some  material  systems  or  applications,  the  process  of  fatigue  crack  formation  in  many  advanced  engineering 
alloys  is  often  complex  and  depends  on  the  influence  of  multiple  interacting  microstructure  attributes.  In  some 
cases,  multiple  interacting  microstructure  attributes  exert  influence  at  different  material  length  scales. 
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Some  recent  computational  work  has  considered  effects  of  multiple  microstructure  attributes  on  fatigue 
response.  Liao  [6]  used  a  Monte  Carlo  technique  to  generate  microstructure  instantiations  with  distributions  of 
particle  size,  grain  size,  and  grain  orientation  that  have  been  randomly  sampled  from  experimental  data  with 
known  distributions.  These  multiple  microstructure  instantiations  are  then  used  to  simulate  the  variation  in  fatigue 
response  in  2024-T351  aluminum  sheets.  The  response  was  estimated  using  FEA  with  elastic  constitutive  laws. 
Fatigue-relevant  attributes,  including  particle  size,  grain  size  and  grain  orientation,  were  identified  in  each 
simulated  instantiation  and  the  distribution  of  these  attributes  was  considered.  These  fatigue-relevant  attributes 
were  selected  a  priori  based  on  experiments.  The  results  of  this  Monte  Carlo  method  correlated  well  with  methods 
based  on  extreme  value  statistics  and  with  the  experiments.  However,  as  with  the  previous  methods  discussed, 
correlations  between  the  important  microstructure  attributes  relative  to  their  influence  on  the  driving  forces  for 
fatigue  crack  formation  were  neglected  because  each  critical  attribute  was  considered  independent  of  the  others. 

1.2  Fatigue  crack  formation  in  duplex  a+p  H-6AI-4V 

In  this  work,  the  processes  of  fatigue  crack  formation  will  be  considered  in  duplex  a+(3  Ti-6AI-4V.  Experimental 
observations  for  Ti-6AI-4V  support  the  hypothesis  that  fatigue  crack  formation  in  this  material  system  is  not 
adequately  described  by  a  distribution  of  any  single  microstructure  attribute.  Fatigue  crack  formation  in  a+(3  Ti 
alloys  has  been  associated  with  slip-dominated  fatigue  crack  formation  as  early  as  the  1960s  by  Wells  and  Sullivan 
[7],  who  observed  cracks  form  along  slip  bands  in  the  a  phase  of  Ti-6AI-4V.  However,  confusion  persists  regarding 
the  fundamental  mechanisms  of  fatigue  crack  formation  in  this  material  system.  In  short,  it  is  not  clear  which 
specific  microstructure  attributes  (or  arrangements  of  attributes)  is/are  most  important  in  the  processes  of  fatigue 
crack  formation.  Relevant  experiments  that  have  investigated  fatigue  crack  formation  in  a+(3  Ti  alloys  under  typical 
cyclical  loading  conditions  at  room  temperature  are  reviewed  next. 

In  duplex  a+P  Ti  alloys,  the  appearance  of  distinct  facets  in  the  primary  a  grains  is  commonly  observed  at  sites 
of  fatigue  crack  formation.  In  some  cases,  these  facets  are  associated  closely  with  basal  planes  oriented 
perpendicular  to  the  uniaxial  loading  axis  in  Ti-6AI-4V  [8,  9].  Bache  [10]  argued  that  crack  formation  on  basal 
planes  in  grains  that  are  unfavorably  oriented  for  slip  is  due  to  shear  stresses  induced  by  dislocation  pileup  on  the 
grain  boundary  in  adjacent  grains  oriented  for  easy  slip.  This  mechanism,  originally  postulated  by  Stroh  [11],  is 
thought  to  activate  when  a  critical  density  of  dislocations  accumulate  on  the  boundary  such  that  the  shear  stresses 
in  the  adjacent  grain  (oriented  for  easy  slip)  are  sufficient  to  induce  slip  on  the  hard  oriented  (i.e.,  nearly 
perpendicular  to  the  loading  direction)  basal  slip  planes.  Thus,  under  cyclic  loading  conditions  a  crack  forms  over 
many  cycles  in  the  hard  oriented  grain.  This  proposed  mechanism  is  supported  by  the  work  of  Baxter  et  al.  [12] 
who  observed  slip  on  basal  planes  oriented  nearly  perpendicular  to  the  loading  axis  in  IMI  834. 

In  contrast,  other  studies  have  observed  fatigue  facet  formation  on  basal  planes  favorably  oriented  for  slip. 
Bridier  et  al.  [13,  14]  investigated  fatigue  crack  formation  in  Ti-6AI-4V  in  both  low  cycle  fatigue  (LCF)  and  HCF, 
respectively,  and  observed  that  cracks  mostly  formed  on  basal  slip  planes  and  less  frequently  on  prismatic  slip 
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planes.  In  all  cases,  these  basal  and  prismatic  planes  exhibited  high  Schmid  factors;  however,  a  bias  was  observed 
towards  primary  a  grains  with  orientation  of  applied  stress  more  towards  the  c-axis,  indicating  a  preference  for 
forming  on  planes  that  are  somewhat  less  favorably  oriented  for  basal  (or  prismatic)  slip  but  with  somewhat  higher 
stress  normal  to  the  basal  (or  prismatic)  plane.  Subsurface  fatigue  crack  formation  on  pyramidal  planes  in  primary 
a  grains  in  Ti-6AI-4V  was  also  observed  by  Gilbert  and  Piehler  [15]  in  grains  oriented  for  hard  basal  slip  (i.e.,  with  c 
axis  nearly  parallel  to  the  uniaxial  loading  direction);  they  noted  that  clusters  of  grains  oriented  favorable  for  slip 
on  pyramidal  <a+c>  planes  and  grains  oriented  favorable  for  basal  slip  tend  to  favor  fatigue  crack  formation.  Jha 
and  Larsen  [16]  observed  in  Ti-6246  that  at  moderate  stress  amplitudes  (i.e.,  close  to  the  macroscopic  yield  stress), 
facets  near  the  location  of  initial  fatigue  crack  formation  were  all  inclined  ~35-45°  relative  to  the  loading  axis.  In 
one  sample,  the  facet  near  the  location  of  initial  fatigue  crack  formation  appeared  to  be  oriented  about  35° 
relative  to  the  loading  direction  in  a  grain  oriented  unfavorably  for  basal  slip,  suggesting  the  facet  may  have  been 
formed  by  a  <c+a>  hard  slip  mechanism  (i.e.,  pyramidal  slip).  In  another  sample,  the  facet  near  the  location  of 
initial  fatigue  crack  formation  developed  on  an  apparent  basal  slip  plane  favorably  oriented  for  slip.  It  is  interesting 
to  note  that  the  grain  unfavorably  oriented  for  basal  slip  that  developed  the  facet  was  surrounded  by  other  grains 
more  favorably  oriented  for  easy  slip  (i.e.,  prismatic  or  basal  slip).  Additionally,  the  sample  with  the  facet  that 
formed  in  the  grain  oriented  for  hard  basal  slip  failed  at  ~60,000  cycles;  the  other  sample  with  the  facet  that 
formed  in  the  grain  oriented  for  easy  basal  slip  failed  at  ~2, 000, 000  cycles.  Szczepanski  et  al.  [17]  also  observed 
fatigue  facet  formation  in  the  very  high  cycle  fatigue  (VHCF)  regime  of  Ti-6246  on  or  near  basal  or  prismatic  slip 
planes  favorably  oriented  for  slip. 

The  size  of  the  primary  a  grains  also  appears  to  play  a  role  in  fatigue  crack  formation  in  duplex  a+P  Ti  alloys. 
Mahajan  and  Margolin  [18]  noted  that  larger  fatigue  cracks  more  readily  formed  in  primary  a  grains  with  longer 
paths  for  slip.  Jha  and  Larsen  [16]  and  also  Szczepanski  et  al.  [17]  noted  in  Ti-6246  that  the  primary  a  grains  that 
form  facets  tend  to  be  larger  than  average.  However,  these  grains  were  not  necessarily  the  largest  primary  a  grains 
in  the  population.  This  suggests  that  grain  size  plays  a  role  in  fatigue  crack  formation  in  this  material  system,  but 
other  attributes  must  also  be  considered. 

The  process  of  fatigue  crack  formation  in  a+P  Ti  alloys  is  further  complicated  when  one  considers  the  existence 
of  local  textured  regions  on  the  scale  of  hundreds  of  microns  to  several  millimeters.  These  microtextured  regions 
are  understood  to  develop  during  processing  from  larger  prior  P  grains.  Bridier  et  al.  [14]  observed  that 
microtextured  regions  in  Ti-6AI-4V  that  favor  slip  activation  in  grains  oriented  for  basal  or  prismatic  slip  are 
associated  with  distributed  fatigue  cracks,  near  the  interface  of  bands  of  texture.  Similar  observations  were  made 
by  Le  Biavant  et  al.  [19],  who  found  distributed  fatigue  cracks  in  microtextured  regions  favorably  oriented  for  basal 
or  prismatic  slip  and  no  cracking  in  regions  unfavorably  oriented  for  prismatic  and  basal  slip.  Because 
microstructure  conditions  in  these  microtextured  regions  appear  to  favor  crack  formation,  multiple  cracks  form  in 
these  regions  that  coalesce  into  a  critical  life  limiting  crack.  Szczepanski  et  al.  [17]  also  noted  that  accumulated  slip 
in  microtextured  regions  with  high  fractions  of  grains  oriented  favorably  for  basal  or  prismatic  slip  in  Ti-6246  likely 
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enhances  fatigue  crack  formation  in  the  VHCF  regime.  Bantounas  et  al.  [20]  observed  that  fatigue  crack  formation 
was  favored  in  regions  with  grains  oriented  for  easy  basal  slip  that  are  adjacent  to  grains  whose  c-axis  are  within 
about  15°  of  the  uniaxial  loading  direction.  They  argued  that  cracks  form  in  grains  that  are  oriented  for  easy  slip 
and  then  propagate  most  effectively  along  basal  planes  in  surrounding  grains  oriented  for  hard  basal  slip.  This  also 
supports  the  previously  mentioned  observations  of  Gilbert  and  Piehler  [15],  who  noted  fatigue  cracking  in  Ti-6AI- 
4 V  in  clusters  of  grains  oriented  for  pyramidal  <a+c>  slip  and  grains  oriented  for  easy  basal  slip. 

Simulations  have  also  given  insight  into  understanding  many  of  these  experimental  observations.  Hasija  et  al. 
[21]  and  Venkataramani  et  al.  [22,  23],  using  dwell  sensitive  crystal  plasticity  models  for  a-Ti-6AI  and  Ti-6246, 
respectively,  demonstrated  that  grains  favorably  oriented  for  basal  slip  transfer  stress  to  nearby  grains  that  are 
unfavorably  oriented  for  basal  slip  through  load  shedding.  Simulations  by  Dunne  [24-26]  for  a-Ti-6AI  showed  that 
these  types  of  rogue  groin  combinations  can  significantly  increase  the  driving  force  for  fatigue  crack  formation, 
particularly  in  dwell  fatigue. 

Based  on  these  previous  examples,  it  is  apparent  that  several  different  microstructure  attributes,  including  grain 
orientation/misorientation,  grain  size  and  local  texture,  have  been  identified  as  important  relative  to  fatigue  crack 
formation  in  a+p  Ti  alloys.  In  fact,  Jha  and  Larsen  [16]  postulate  a  hierarchy  of  microstructure  arrangements 
relative  to  their  potency  for  localized  accumulated  plastic  slip.  Based  on  limited  experiments,  they  suggest  four 
possible  underlying  microstructure  configurations  in  Ti-6246  whose  probability  of  occurrence  may  relate  to  the 
observed  variability  in  the  overall  fatigue  life.  A  cluster  of  primary  a  grains  and  lamellar  a+P  colonies  with  aligned 
slip  planes  oriented  for  close  to  maximum  shear  or  a  cluster  of  primary  a  grains,  each  similarly  oriented  for  basal 

(aj  slip  were  observed  at  sites  of  subsurface  fatigue  crack  formation  in  specimens  with  shorter  fatigue  lives.  In 

contrast,  the  microstructure  features  described  by  a  hard  primary  a  grain  (oriented  unfavorably  for  basal  ^ a ^  slip) 
surrounded  by  soft  lamellar  colonies  and/or  primary  a  grains  favorably  oriented  for  slip  or  primary  a  grains 
favorably  oriented  for  basal  ^ a ^  slip  were  observed  to  be  associated  with  surface  fatigue  crack  formation  in 

specimens  with  longer  fatigue  lives.  Likely,  many  more  experiments  would  be  needed  to  make  any  conclusions 
regarding  the  uniqueness  of  this  set  of  critical  microstructure  configurations  for  fatigue  damage  formation  in  this 
material  system. 

This  body  of  work  indicates  that  the  mechanisms  of  fatigue  crack  formation  in  these  types  of  a+P  Ti  alloys  are 
quite  complex.  Even  in  HCF,  localized  plastic  slip  is  clearly  linked  to  fatigue  crack  formation.  Additionally,  the 
multiplicity  of  competing  failure  modes  increase  fatigue  variability,  particularly  in  the  HCF  and  VHCF  regimes  [27- 
29].  It  is  not  likely  that  the  single  attribute  extreme  value  approaches  outlined  previously  will  be  able  to  capture 
the  microstructure  dependence  of  the  fatigue  crack  formation  and  growth  in  these  types  of  material  systems.  A 
new  approach  is  needed  to  link  extreme  value  distributions  of  interacting  microstructure  attributes  that  enables 
the  identification  and  comparison  of  configurations  of  microstructure  attributes  that  drive  variability  in  fatigue  life. 
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1.3  Correlated  microstructure  attributes  and  extreme  value  driving  forces  for  fatigue  crack  formation 

Recently  a  new  probabilistic  framework  was  introduced  that  uses  extreme  value  marked  correlation  functions 
to  quantify  the  interacting  microstructure  attributes  that  are  most  statistically  significant  relative  to  the  measured 
or  simulated  extreme  value  distributions  of  response  [30,  31].  Specifically,  correlation  functions  describing  the 
microstructure  attributes  potentially  relevant  to  the  operant  mechanism  of  fatigue  crack  formation  are  selectively 
sampled  in  regions  marked  by  certain  extreme  value  response  parameters  (e.g.  cyclic  plastic  strain  range).  These 
extreme  value  marked  correlations  can  be  compared  to  the  same  correlation  functions  sampled  over  the  entire 
microstructure  ensemble  (without  regard  to  response)  to  determine  the  particular  correlated  microstructure 
attributes  that  are  most  significant  to  the  extreme  value  response.  In  general,  this  methodology  can  be  applied  to 
quantify  important  microstructure  attributes  that  correlate  with  strong  influence  on  the  extreme  value 
distributions  of  fatigue  response.  This  probabilistic  framework  was  coupled  with  multiple  simulated  material 
volumes  to  investigate  orientation  relationships  important  to  fatigue  crack  formation  in  the  PM  Ni-base  superalloy 
IN  100  [30,  31].  It  was  shown  that  the  driving  forces  for  fatigue  crack  formation  were  highest  in  clusters  of  grains 
oriented  for  cube  slip  near  grains  oriented  for  octahedral  slip,  which  correlated  well  with  experimental 
observations  in  the  same  material  [32]. 

Simulation-based  strategies  for  investigating  the  driving  forces  for  fatigue  crack  formation  at  the  scale  of 
microstructure  have  greatly  expanded  in  recent  years.  McDowell  [33]  outlined  the  use  of  plastic  shear  strain- 
based  FIPs  to  estimate  the  driving  forces  for  fatigue  crack  formation  in  metals  associated  with  plastic  ratcheting 
and  reversed  cyclic  plasticity,  i.e., 


end  of  the  cycle 


beginning  of  the  cycle 


(i) 


and 


end  of  the  cycle 


(2) 


respectively,  where  S?  are  components  of  the  plastic  strain  tensor.  Figure  1  illustrates  these  cases.  Two  such 
parameters  are  the  maximum  plastic  shear  strain  range  and  Fatemi-Socie  [34]  damage  parameter,  i.e., 
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The  nonlocal  maximum  cyclic  plastic  shear  strain  range  is  computed  using  the  ordered  principal  plastic  strain 
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employed  to  denote  averaging  of  a  quantity  "  *  "  over  a  characteristic  volume,  V.  Note  that  averaging  is  performed 
on  each  component  of  the  plastic  strain  tensor,  i.e.  (zfj^J  =  (l/K)  J  sfjdV  •  The  peak  stress  normal  to  the  plane 

of  Xmax)K  's  S'ven  hy  cr^ax,  cry  is  macroscopic  cyclic  yield  strength,  and  K  controls  the  effect  of  the  normal 

stress  on  fatigue  crack  formation.  Here  a  value  of  K  =  l  is  used  based  on  previous  studies  (e.g.,  [31]),  and  is 
subject  to  refinement.  Averaged  over  scales  pertinent  to  the  microstructure  attributes  relevant  to  the  processes 
of  fatigue  crack  formation,  these  types  of  FIPs  have  been  used  to  investigate  the  relative  fatigue  crack  formation 
potency  of  several  different  microstructure  attributes.  In  martensitic  gear  steels,  Prasannavenkatesan  et  al.  [35] 
used  these  two  parameters  to  characterize  the  relative  potency  of  different  kinds  of  inclusions  to  form 
grain/inclusion  scale  HCF  cracks.  Similar  FIPs  have  also  been  used  to  characterize  the  influence  of  microstructure 
on  fatigue  crack  formation  in  structural  steel  [36,  37]  and  Ni-base  superalloys  [25,  30,  31,  38,  39].  These  efforts 
have  demonstrated  the  utility  in  using  these  types  of  parameters  coupled  with  crystal  plasticity  simulations  to 
correlate  microstructure  scale  slip  with  fatigue  crack  formation  at  the  scale  of  individual  grains  or  attributes  that 
control  slip  intensification  (e.g.,  inclusions).  Although  not  pursuant  to  the  present  work,  It  is  possible  to  confer 

grain  scale-specific  Coffin-Man  son  type  laws  for  crack  formation,  e.g.,  PFS  =  y'f  ( Ninc  )c ,  where  is  a  rescaled 

fatigue  ductility  coefficient  and  C  is  the  fatigue  ductility  exponent. 

In  this  work,  the  same  microstructure-sensitive  extreme  value  probabilistic  framework  is  used  that  was 
demonstrated  previously  for  IN  100  [30,  31]  to  investigate  the  influence  of  the  microstructure  on  the  driving  forces 
for  fatigue  crack  formation  in  duplex  Ti-6AI-4V.  A  similar  strategy  of  simulating  multiple  statistical  volume  elements 
(SVEs)  sampled  from  the  same  distributions  of  microstructure  is  used  to  characterize  the  extreme  value 
distributions  of  FIPs.  Extreme  value  marked  correlation  functions  are  sampled  selectively  at  the  locations  of 
extreme  value  response  (i.e.,  hot  spots).  These  statistical  volume  elements  are  simulated  via  the  finite  element 
method  using  a  rate  sensitive  crystal  plasticity  model  calibrated  to  experimental  measurements. 

2.  Methodology 

2.1  Extreme  value  marked  microstructure  correlation  function 

A  brief  review  of  the  previously  developed  microstructure-sensitive  extreme  value  probabilistic  framework  is 
given  [30,  31].  Based  on  classical  extreme  value  probability  theory  (cf.  [4]),  the  probability  distribution  of  the 
extreme  valued  response  parameter  a  (e.g.,  a  =  FIP)  is  defined  as  Fex  (a  |Q)  given  a  window  (i.e.,  SVE)  Q. 
Specifically,  Fex  ( a  |  Q)  describes  the  probability  that  the  response  parameter  of  value  a  is  the  extreme  value  for  a 

sampled  SVE  of  microstructure  Q .  The  extreme  value  marked  radial  correlation  function  is  then  defined  as 
Rmax(a)  |  r,Q)dr ,  denoting  the  probability  that  a  microstructure  attribute  /3  located  coincident  with  the 

extreme  value  response  parameter  a  in  Q  is  located  within  a  distance  of  r  to  r  +  dr  of  a  second  microstructure 
attribute  /?'  in  any  direction.  An  extreme  value  marked  n-point  correlation  function,  extreme  value  marked  lineal 
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path  correlation  function,  extreme  value  marked  nearest  neighbor  correlation  function,  etc.  can  similarly  be 
defined  [40].  This  two  part  statistical  construct  describes  both  the  extreme  value  response  of  the  microstructure  as 
represented  by  the  response  parameter  OC  and  considers  the  probabilities  of  correlated  microstructure  attributes, 
1 3  and  /?',  existing  in  the  proximity  of  extreme  values  of  a  in  a  microstructure  window  Q.  By  comparing 
R  max  («)  ^  ^  |  r?Q)jr  with  the  correlation  function  for  the  same  microstructure  attributes  sampled  from  the 
complete,  bulk  material  ensemble,  given  by  R  (/?,/?'  |  r)  dr ,  the  correlated  microstructure  attributes  most 
probable  to  exist  at  locations  of  extreme  value  response  can  be  identified.  Sampling  is  performed  on  a  predefined 
size  of  microstructure  window  Q  <  Changing  the  size  of  the  SVE  affects  the  number  of  samples  required  to  build 
up  the  statistics  of  the  tail  of  the  extreme  value  distribution;  fewer  are  required  as  SVE  size  increases.  Moreover, 
the  minimum  size  of  Q  is  established  by  relevant  spatial  correlations  of  the  neighborhood  of  extreme  value 
response  sites.  For  example,  nearest  and  second  nearest  neighbor  grains  must  be  sampled  to  capture  load 
shedding  phenomena,  while  much  larger  values  are  perhaps  required  to  capture  phenomena  associated  with 
microtexture.  In  other  words,  the  SVE  size  must  be  suitably  long  compared  to  the  short  range  spatial  interactions 
at  microstructure  hot  spots.  In  general,  if  the  SVE  size  is  too  small,  the  extreme  value  distribution  resulting  from 
compiling  extreme  values  of  a  set  of  SVEs  may  not  relate  to  the  low  probability  tail  region  of  interest,  and  may 
provide  ambiguous  information. 

2.2  Crystal  Plasticity  Modeling  ofTi-6Al-4V 

In  aerospace  applications,  a+P  Ti  alloys  are  important  in  view  of  their  high  strength  and  desirable  strength-to- 
weight  ratio.  These  alloys  exhibit  favorable  properties  at  elevated  (but  moderate)  temperatures  depending  on  heat 
treatment  and  thermomechanical  processing.  The  microstructure  of  typical  a+P  Ti  alloys  can  range  from  bi-modal 
to  fully  lamellar.  A  bi-modal  microstructure  is  generated  by  applying  deformation  when  the  material  is  in  the  a+P 
regime  (on  a  phase  diagram)  followed  by  recrystallization  and  ageing.  This  processing  route  generates  a 
microstructure  that  consists  of  equiaxed  primary  a  grains  dispersed  among  regions  containing  transformed  (3 
grains  with  embedded  a  laths.  Texture  effects  in  these  alloys  are  particularly  pronounced  due  to  the  low  symmetry 
of  crystallographic  slip  and  anisotropic  elasticity,  as  well  as  the  contrast  of  properties  between  the  a  and  p  phases. 
In  addition,  despite  their  greater  number  of  available  slip  systems,  the  P  regions  containing  a  laths  tend  to  be  more 
resistant  to  slip  than  the  primary  a  grains  due  to  the  presence  of  the  a  laths  and,  in  some  cases,  finer  a 
precipitates. 

The  particular  crystal  plasticity  model  used  in  this  work  for  duplex  Ti-6AI-4V  was  initially  developed  by  Mayeur 
and  McDowell  [41],  extended  by  Zhang  and  McDowell  [42]  and  then  Bridier  et  al.  [43]  specifically  for  application  to 
the  HCF  loading  regime  for  peak  stresses  at  or  below  the  macroscopic  yield  strength.  In  hep  primary  a,  there  are 
four  different  families  of  slip  systems:  three  ^1120^(0001)  basal,  three  ^1120)(l010|  prismatic,  six 
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^1120)jl01l}  first  order  pyramidal  <a>,  and  twelve  ^1123)jl01l}  first  order  pyramidal  <a+c>  slip  systems. 

However,  slip  is  dominant  in  the  basal  and  prismatic  slip  systems  due  to  a  relatively  low  critical  resolved  shear 
stress  compared  to  the  other  slip  systems.  The  lamellar  a+(3  colonies  are  homogenized  in  this  model  since  the 
individual  laths  range  in  thickness  from  hundreds  of  nanometers  to  several  microns.  A  crystallographic  burgers 
orientation  relation  (BOR)  is  maintained  between  the  secondary  a  and  (3  laths  defined  such  that  (000l)a//  {101}. 

and  (000l)  //(lll)p.  The  24  possible  slip  systems  in  the  lamellar  region  include  three  ^1120^(0001)  basal,  three 

^1120)(l010|  prismatic,  six  ^1120)(l01l|  first-order  pyramidal  and  twelve  (lll){ll0}  bcc  slip  systems.  The  bcc 

slip  systems  are  transformed  into  the  hexagonal  coordinate  system  via  the  BOR.  Hard  slip  systems  in  these  colonies 
are  those  that  intersect  the  a/(3  interface.  Soft  deformation  modes  are  those  on  which  dislocations  glide  parallel  to 
the  a/P  interface  or  have  parallel  slip  planes  in  both  the  secondary  a  and  P  phases.  In  particular,  Bridier's  [43] 
model  is  formulated  to  favor  single  slip,  which  has  been  experimentally  observed  to  dominate  at  low  cyclic  strain 
amplitudes.  In  addition,  slip  in  the  a/P  colonies  is  more  restricted. 

2.3  Ellipsoidal  Packing  Algorithm  for  Generation  of  Digital  SVEs  for  Ti-6Al-4V 

The  extreme  value  fatigue  response  in  Ti-6AI-4V  is  characterized  via  the  simulation  of  multiple  SVEs  to  estimate 
Fex(a  |Q)  and  Rmtt*W  (p9p’ \r,Cl)-  Cyclic  loading  response  is  computed  using  3D  digital  SVEs  of  microstructure 

using  the  FE  package  ABAQUS  [44]  coupled  with  a  crystal  plasticity  model  calibrated  with  experiments  on  duplex 
Ti-6AI-4V,  implemented  in  a  user-material  subroutine.  Using  multiple  SVEs  sampled  from  predefined  parent 
distributions  of  key  microstructure  attributes  (e.g.  grain  size  distribution,  phase  distribution,  orientation 
distribution,  misorientation/disorientation  distribution,  etc.),  along  with  3D  periodic  boundary  conditions,  allows 
the  simulation  of  microstructure  dependent  variation  in  fatigue  response. 

The  size  of  the  SVE  for  a  desired  response  (i.e.,  grain  scale  plasticity)  is  selected  to  be  large  enough  such  that  the 
response  parameter  of  interest  at  a  particular  location  is  unaffected  by  statistical  variations  in  the  microstructure 
at  a  distance  away  on  the  order  of  the  SVE  size;  however,  the  volume  need  not  be  so  large  that  it  contains  the 
entire  distribution  of  response  for  that  particular  response  parameter.  This  is  in  contrast  to  a  representative 
volume  element  (RVE),  defined  to  be  sufficiently  large  to  ensure  that  the  distribution  of  local  response  or  effective 
response  of  the  volume  will  not  change  by  increasing  the  size  of  the  RVE  or  its  position  in  the  bulk  material. 
However,  a  RVE  of  reasonable  size  is  often  impractical,  either  being  larger  than  laboratory  specimens  or  involving 
far  too  intensive  calculations.  This  is  particularly  true  when  trying  to  capture  extreme  value  distributions  of 
microstructure  response.  In  this  work,  the  size  of  the  SVEs  is  chosen  such  that  the  local  plastic  deformation 
response  on  the  scale  of  the  grains  is  not  substantially  affected  by  further  increasing  SVE  size.  In  addition  to  local 
variations  in  slip  response,  it  is  also  desired  that  these  volumes  be  representative  (i.e.,  be  an  RVE)  in  terms  of  the 
overall  effective  (i.e.,  macroscopic)  elastic/plastic  response.  In  other  words,  it  is  desired  that  the  stress-strain 
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response  over  the  SVE  taken  from  any  microstructure  instantiation  to  be  representative  of  any  other.  This  is  the 
same  definition  of  statistical  volume  element  used  in  previous  work  [31]  and  very  similar  to  the  definition  of 
SERVES  introduced  by  Swaminathan  and  Ghosh  [45,  46]. 

To  generate  SVEs  Ti-6AI-4V,  an  ellipsoidal  packing  based  algorithm  is  employed  to  generate  digital 
microstructure  models.  This  method  uses  known  distributions  of  the  aspect  ratios  of  grain  equivalent  ellipsoids  fit 
to  grains  of  experimentally  characterized  microstructures.  Using  grain  equivalent  ellipsoids  instead  of  other  space 
filling  methods  like  Voronoi  tessellation  allows  construction  of  more  complex  grain  morphologies  such  as 
elongated  grains  common  in  rolled  ductile  metals,  and  facilitates  more  accurate  reconstruction  of  two  phase 
microstructures  with  bi-modal  size  distribution  such  as  duplex  Ti-6AI-4V.  For  example,  Brahme  et  al.  [47]  used 
ellipsoids  fit  to  actual  grain  distributions  to  reconstruct  Al  microstructures  with  elongated  grains.  The  algorithm  use 
here  was  developed  in  part  based  on  the  work  of  Groeber  et  al.  who  developed  ellipsoidal  based  methods  to 
reconstruct  polycrystalline  representations  of  IN  100  characterized  from  three-dimensional  serial  sections  [48,  49]. 

The  ellipsoidal  packing  based  digital  microstructure  generation  algorithm  was  implemented  in  MATLAB  [50], 
with  logic  as  follows: 

1)  An  ellipsoid  is  generated  by  randomly  selecting  the  ellipsoid  aspect  ratios  from  a  predefined  distribution 
(e.g.,  beta  distribution). 

2)  The  orientation  of  the  ellipse  is  randomly  selected  based  on  a  predefined  distribution  (i.e.,  orientation 
distribution  function  for  orientation  of  major  axes  of  grain  equivalent  ellipsoids  relative  to  the  specimen 
axis). 

3)  The  ellipsoid  is  randomly  placed  in  the  volume.  If  the  ellipsoid  overlaps  with  any  previously  placed 
ellipsoid,  a  new  random  location  is  selected  repeatedly  until  no  overlap  exists  (see 

4)  Figure  2a. 

5)  Steps  1  through  4  are  repeated  until  the  jamming  limit  is  reached  and  no  more  non-overlapping  ellipsoids 
can  be  placed  in  the  volume. 

6)  The  ellipsoids  are  then  allowed  to  grow  uniformly  until  all  the  space  in  the  microstructure  volume  is  filled. 

A  two-dimensional  section  of  a  generated  microstructure  using  this  method  with  a  bi-modal  grain  size 
distribution  is  shown  in  Figure  2b.  Comparing  the  aspect  ratios  of  the  ellipsoids  and  final  grains  after  applying 
uniform  grain  growth,  the  distributions  of  the  aspect  ratios  of  the  original  ellipsoids  and  the  equivalent  ellipsoids  fit 
to  the  final  grains  were  essentially  the  same,  as  shown  in  Figure  3.  This  method  gives  the  user  much  more  control 
over  grain  morphology  than  other  commonly  used  microstructure  generation  algorithms  like  those  based  on  the 
Voronoi  tessellation  and  allows  one  to  easily  generate  bi-modal  grain  size  distributions  as  shown  in  Figures  4  and  5. 
Here,  the  tails  of  the  distributions  are  captured  with  higher  fidelity  than  the  tails  of  the  grain  size  distribution  for 
the  microstructures  generated  for  IN  100  previously  using  a  Voronoi  tessellation  based  microstructure  generator 
[31]. 
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2.4  Simulated  Microstructures 

To  explore  the  effects  of  varying  microstructure  on  the  distributions  of  FIPs,  four  microstructures  were 
generated  with  different  grain  size  distributions  and  volume  fractions  of  primary  a.  The  target  size  of  the  primary  a 
grains,  a+(3  colonies,  and  the  target  volume  fraction  of  the  primary  a  phase  are  given  in  Table  1.  Grain  size  is 
related  to  grain  volume  according  to 

K)-0.1{vT  (5) 

where  (dgrT^j  and  (V^  are  the  mean  (volume  averaged)  grain  size  and  grain  volume,  respectively  [51].  The 

assumed  mean  (jli)  and  standard  deviation  (a)  of  the  grain  size  for  these  microstructures  is  also  given  in  Table  1.  In 
Figures  4  and  5,  the  target  log-normal  and  fit  cumulative  distribution  functions  are  shown  for  the  primary  a  grains 
and  a+P  colonies  for  randomly  generated  fine  and  coarse  bi-modal  microstructures,  respectively. 

2.3  Extreme  Value  Statistics  of  Fatigue  in  H-6AI-4V 

This  work  considers  how  the  crystallographic  attributes  of  the  microstructure  (e.g.,  grain  orientation,  grain 
disorientation,  grain  size,  and  grain  shape  distributions)  affect  local  driving  forces  for  fatigue  crack  formation  in 
HCF.  As  introduced  previously,  the  maximum  plastic  shear  strain  (MPSS)  and  Fatemi-Socie  (FS)  parameters  are 

employed.  Additionally,  the  cumulative  effective  plastic  strain  (CEPS),  s?ff  =  yj( 2/3)^^  ,  over  all  cycles  is 

considered.  In  each  microstructure,  these  parameters  are  calculated  using  each  of  three  different  sized  averaging 
volumes,  as  defined  in  Table  3.  The  smallest  averaging  volume  considered  is  slightly  larger  than  the  mean  grain  size 
of  the  smallest  grains  (i.e.,  the  primary  a  grains  in  Microstructures  A  and  B).  By  using  the  same  sized  averaging 
volumes  to  compare  the  response  across  the  range  of  simulated  microstructure,  a  uniform  comparison  can  be 
made  of  the  driving  forces  for  fatigue  crack  formation  on  the  scale  of  the  averaging  volume  regardless  of  the 
differing  distributions  of  grains/phases  sizes. 

Multiple  SVEs  are  constructed  for  Ti-6AI-4V  using  the  ellipsoidal  based  microstructure  generator  and  then 
simulated  via  FEM.  In  all  cases,  these  simulations  are  loaded  with  a  positive  mean  stress.  These  SVEs  are  all 
dimensioned  0.400  mm  along  each  edge  of  the  cuboidal  microstructure  block  and  are  subjected  to  uniaxial  cyclic 
strain  at  0.6%  peak  strain  with  zero  minimum  strain  (R8  =  0)  at  a  quasistatic  strain  rate  of  0.002s  1  at  room 
temperature.  Each  simulation  is  cycled  at  least  ten  and  up  to  twenty  times  to  explore  how  shakedown  affects  the 
extreme  value  distributions  of  the  FIPs.  Periodic  boundary  conditions  are  applied  in  all  directions  to  simulate 
subsurface  conditions.  The  voxellated  meshes  of  these  generated  microstructures  consist  of  quadrilateral  elements 
with  reduced  integration  (type  C3D8R  in  ABAQUS).  The  specific  range  of  the  number  of  grains  is  given  in 

Table  2  for  the  different  microstructures  considered.  There  are  26  elements  along  each  edge  of  the  SVE  or  17576 
elements  in  all.  The  volume  of  the  SVE  was  maintained  constant  in  this  analysis  across  the  range  of 
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microstructures.  A  meshed  SVE  for  an  arbitrary  instantiation  of  Microstructure  A  is  given  in  Figure  6  with  contours 
of  the  Mises  strain  at  the  peak  strain  Of  0.6%  after  10  cycles.  The  variations  in  the  stress  fields  shown  in  Figure  6 
throughout  the  volume  are  primarily  due  to  the  mismatch  of  the  anisotropic  elastic  stiffness  for  the  differently 
oriented  grains. 


3.  Results  and  Discussion 

To  understand  the  number  of  cycles  required  to  approach  a  condition  of  plastic  shakedown  (i.e.,  stabilization  of 
cyclic  plasticity  at  the  grain  scale),  a  single  instantiation  for  Microstructure  A  was  cycled  20  times  using  the 
conditions  given  previously.  The  percentage  of  primary  a  grains  and  a+(3  colonies  that  exhibit  an  effective 
cumulative  plastic  strain  above  a)  l.OxlO’10,  b)  l.OxlO’8,  and  c)  l.OxlO’6  are  shown  in  Figure  7.  Based  on  the 
shakedown  period  in  Figure  7,  it  was  assumed  that  at  least  10  cycles  were  required  before  relative  stability  was 
achieved  for  the  cyclic  plastic  response.  Although  more  cycles  are  desirable,  due  to  computational  limitations  all 
subsequent  simulations  were  subjected  to  10  cycles  to  illustrate  application  of  the  algorithms. 

Next,  the  extreme  value  distributions  of  the  various  FIPs  described  above  are  considered  for  the  four  different 
microstructures  described  in  Table  1.  A  total  of  100  SVEs  were  instantiated  for  each  microstructure  in  Table  1.  The 
FIPs  were  calculated  over  three  different  averaging  volumes  as  defined  in  Table  3  along  with  the  equivalent  grain 
size  for  those  averaging  volumes  as  estimated  via  Equation  (5).  A  volume  average  FIP  centered  at  each  element  in 
each  SVE  was  computed  for  each  averaging  volume.  The  maximum  volume  averaged  FIP  for  a  given  SVE  was 
defined  as  the  extreme  value  FIP  for  that  instantiation.  The  distribution  of  extreme  value  FIPs  was  then 
constructed  for  each  microstructure  based  on  values  from  each  of  100  SVEs.  For  the  smallest  averaging  volume, 
these  distributions  are  plotted  in  Figure  8  on  a  Gumbel  probability  scale  [31].  The  distribution  appears  linear  on 
the  Gumbel  scale  if  it  fits  that  particular  distribution  exactly.  Using  linear  least  squares  regression,  the  distribution 
parameters  for  the  particular  distribution  of  interest  can  be  estimated. 

Comparing  the  extreme  value  distributions  for  three  different  FIPs  does  not  immediately  reveal  any  significant 
differences  between  the  distributions;  the  same  general  relative  trends  are  observed  across  the  range  of 
microstructures.  However,  the  grains  predicted  to  exhibit  the  extreme  value  response  varied  between  the 
different  FIPs  in  as  many  as  10%  of  the  SVEs  simulated.  It  is  likely  that  CEPS  is  the  least  desirable  FIP,  as  the 
cumulative  effective  plastic  strain  can  evolve  dramatically  during  initial  cycles,  but  then  changes  little  during 
subsequent  cycling.  Unless  it  can  be  definitively  shown  that  directional  plastic  strain  accumulation  drives  fatigue 
crack  formation  in  very  early  cycles  (unlikely  in  HCF),  a  measure  of  steady  state  cycle  by  cycle  plastic  strain  reversal 
or  accumulation  is  preferred.  The  FS  FIP  likely  best  captures  the  actual  mechanisms  of  fatigue  crack  formation 
based  on  the  previous  experimental  observations  of  the  local  damage  processes  acting  in  this  class  of  materials, 
although  CEPS  has  been  used  in  the  literature  [24-26].  The  FS  parameter  is  better  suited  to  capture  the  effect  of 
the  normal  stress  that  appears  to  promote  formation  of  fatigue  cracks  on  basal  or  prismatic  slip  planes  with  slip 
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plane  normals  that  tend  to  be  closer  aligned  with  the  loading  direction  [43].  As  shown  by  McDowell  and  Berard 
[52],  the  FS  parameter  reflects  mixed  mode  growth  of  small  fatigue  cracks.  For  Ti-6AI-4V  it  is  expected  to  be 
effective  in  reflecting  combined  effects  of  basal  slip  with  elevated  normal  stress  to  this  plane  due  to  intergranular 
interactions. 

As  discussed  in  previous  work  [31],  it  is  of  interest  to  consider  the  mathematical  character  of  the  extreme  value 
distributions  of  the  FIPs  as  estimated  from  the  simulated  SVEs.  Assuming  that  X  is  a  random  variable,  the 
maximum  extreme  value  of  a  sample  from  X  of  size  n  can  be  defined  as 

Yn=max(X1,X2,...,Xn).  (6) 

If  it  is  required  that  Yn  is  less  than  some  value  y ,  then  all  the  random  variables  in  the  same  sample  associated 
with  Yn  must  also  be  less  than  y .  The  distribution  function  of  Yn  is  defined  as  the  probability  that  for  a  sample  of 


size  n  ,  Yn  is  less  than  or  equal  to  y  ,  i.e., 


Fyn  (. y )  -  P(Y„  <  y)  =  P(X,  <  <  y,...,Xn  <  y ) • 


(7) 


For  distributions  of  a  single  variable,  it  has  been  shown  that  there  are  only  three  types  of  non-degenerated 
distributions  to  which  the  extreme  value  distributions  can  converge  for  sufficiently  large  sample  sizes  [53,  54].  The 
three  possible  non-degenerated  asymptotic  distributions  for  the  maximum  extreme  value  distributions  are  the  (i) 
Gumbel  (Type  I),  (ii)  Frechet  (Type  II),  or  (ii)  Weibull  (Type  III).  Mathematically,  these  distributions  can  take  the 
respective  forms  [55] 


FYH(y„)=QxP 


(8) 


fy„  (y„ )  =  exP 


and 


FY„  (>’„)  =  exP 


(  Yl 

UJ  _ 

f 

c°~y„ 

{®-wn 

(9) 


(10) 


For  the  Gumbel  distribution,  un  is  the  characteristic  largest  value  of  the  initial  variable  X  and  an  is  an  inverse 

measure  of  dispersion  of  the  largest  value  of  X .  The  number  of  samples  of  the  initial  variable  X  in  the  set  of 
distributions  of  X  from  which  Yn  is  sampled  is  n  .  For  the  Frechet  distribution,  vn  is  the  characteristic  largest 


value  of  the  underlying  variable  X  and  k  is  a  shape  parameter  or  measure  of  the  dispersion.  Finally,  in  the 
Weibull  distribution,  co  is  the  upper  bound  of  the  initial  parent  distribution  (i.e.,  Fx^(o)  =  1),  wn  is  the 

characteristic  largest  value  and  k  is  the  shape  parameter.  The  desire  is  to  characterize  a  dataset  as  belonging  to 
one  of  these  three  types  of  distributions  so  that  mathematical  properties  can  be  used  to  better  understand  the 
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data.  Here  the  Type  III  Weibull  distribution  is  not  considered  because  the  authors  at  this  time  are  unaware  of  a 
method  to  quantifying  an  upper  bound  co  for  the  present  data. 

In  Figure  9,  the  extreme  value  FS  FIPs  calculated  over  a  cube  shaped  averaging  volume  with  equivalent  grain  size 
of  0.032mm  are  plotted  on  the  a)  Gumbel  (Type  I)  and  b)  Frechet  (Type  II)  probability  scale  along  with  the  linear 
fits  estimated  via  least  squares  regression.  The  Frechet  probability  scale  is  constructed  in  the  same  manner  as  the 
Gumbel  probability  scale  that  is  described  in  [31].  The  linearized  Frechet  distribution  is  of  the  form 

ln[l/ln(l/p)]  =  k\n(yn)-k\n(vn)  ,  where  p  is  the  probability.  The  fitting  parameters  for  the  maximum 
extreme  value  Gumbel  distribution  defined  by  Equation  (8)  and  the  Frechet  distribution  defined  by  Equation  (9) 
are  given  in  Tables  4  and  5.  According  to  the  calculated  R 2  values  from  the  least  squares  regression  of  the  entire 
extreme  value  dataset,  the  distributions  appear  to  best  fit  the  Frechet  distribution.  However,  upon  visual 
inspection,  it  is  obvious  that  the  distributions  are  not  linear  in  Figure  9b.  In  fact,  the  distributions  are  much  more 
linear  on  the  Gumbel  probability  scale  in  Figure  9a,  except  for  the  lower  tails  (below  about  0.5  probability)  which 
tend  to  have  a  steeper  slope  then  rest  of  the  distribution.  The  higher  tails  (0.99  probability  and  above)  for 
Microstructure  B  and  D  also  exhibit  increased  scatter  and  deviate  significantly  from  the  more  linear  portion  of  the 
distributions  on  the  Gumbel  probability  scale.  Here  there  appears  to  be  rare  cases  when  the  volume  fraction  of  the 
primary  a  grain  is  much  higher  (i.e.,  70%  as  in  the  case  of  Microstructure  B  and  D)  that  there  are  extreme  value 
FIPs  that  appear  to  jump  up  significantly  magnitude  relative  to  the  rest  of  the  distribution.  It  could  be  the  case  that 
different  regimes  of  the  FIPs  associate  with  lower  and  higher  probabilities  of  the  extreme  value  dataset  In  terms  of 
microstructure  interactions;  this  bears  further  examination  in  future  work.  However,  in  the  present  work  it  is  the 
Gumbel  distribution  appears  to  best  characterize  most  of  the  distribution  with  the  noted  exceptions. 

In  Figure  10,  the  extreme  value  distribution  of  the  FS  FIP  determined  for  each  of  the  first  10  cycles  for 
Microstructure  A  is  given  for  the  FIPs  calculated  over  a  cube  shaped  averaging  volume  with  equivalent  grain  size  of 
0.032mm.  From  Figure  10,  it  can  be  observed  that  most  of  the  shakedown  occurs  during  the  first  6-7  cycles;  there 
is  little  change  in  the  distributions  for  cycles  8-10.  Shakedown  of  the  extreme  value  distributions  can  be  quantified 
by  comparing  the  parameters  of  the  linear  fits  for  the  Gumbel  distribution  that  are  extracted  as  previously 
explained.  The  inverse  measure  of  the  dispersion  of  these  distributions,  an,  changes  by  as  much  as  50%  between 

the  initial  and  second  cycles  but  changes  by  less  than  6%  between  the  cycles  9  and  10.  However,  the  grains 
predicted  to  exhibit  the  extreme  value  FS  FIPs  differed  less  than  3%  of  the  time  between  the  1st  and  10th  cycles, 
indicating  that  location  of  the  extreme  value  FIP  in  most  cases  can  be  detected  in  early  cycles. 

In  Figure  11,  the  distributions  of  the  volume  average  Fatemi-Socie  (FS)  FIP  calculated  over  an  averaging  volume 
with  equivalent  grain  size  of  a)  0.054mm  and  b)  0.076mm  are  shown  on  a  Gumbel  probability  scale.  Comparing 
Figures  9a  and  11,  there  is  no  significant  difference  in  the  relative  distributions  of  the  extreme  value  FS  FIPs 
between  the  different  averaging  volumes  considered.  The  larger  the  averaging  volume  for  the  FIP  calculations,  the 
lower  the  overall  magnitude  of  the  extreme  value  FIPs  tends  to  be.  This  is  reasonable  because  larger  averaging 
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volumes  will  wash  out  the  extreme  value  response  over  a  larger  area.  Using  larger  averaging  volumes,  however, 
allows  one  to  consider  the  average  driving  forces  for  fatigue  crack  formation  over  a  larger  area.  Depending  on  the 
local  barriers  for  microstructurally  small  crack  growth,  larger  averaging  volumes  (relative  to  the  distance  between 
the  dominant  microstructural  barriers)  might  better  indicate  the  locations  where  the  critical  life  limiting  cracks 
form  that  are  large  enough  to  avoid  being  arrested  by  these  barriers. 

Based  on  the  extreme  value  distributions  of  the  FIPs,  Microstructure  A  is  predicted  to  have  the  least  variability 
in  overall  fatigue  performance  and  lower  maximum  driving  forces  for  fatigue  crack  formation.  Thus,  it  is  likely  that 
microstructures  with  small  primary  a  grain  sizes  and  smaller  volume  fractions  of  the  primary  a  grains  will  have 
better  fatigue  performance.  This  is  directly  related  to  the  fact  that  most  cyclic  slip  is  occurring  in  the  primary  a 
grains  at  the  applied  peak  strain  of  0.6%.  Lower  volume  fractions  of  the  primary  a  grains  exhibit  less  slip  activity 
than  the  microstructures  with  higher  volume  fractions  of  primary  a  grains.  In  addition,  the  smaller  primary  a  grains 
have  smaller  mean  free  paths  for  slip. 

The  microstructure  attributes  that  are  coincident  with  the  grains  in  which  the  extreme  value  FIPs  were 
identified  are  considered  next.  All  of  the  volume  averaged  extreme  value  FIPs  calculated  over  a  cube  shaped 
averaging  volume  with  equivalent  grain  size  of  0.032mm  were  identified  to  be  centered  within  primary  a  grains.  To 
consider  the  orientation  of  the  grains  predicted  to  exhibit  the  estimated  extreme  value  response,  Schmid  factors 

were  identified  for  the  primary  a  HCP  slip  systems  including  the  basal,  prismatic,  pyramidal  ,  and  pyramidal 

(a  +  cj  slip  systems.  Figure  12  shows  the  apparent  Schmid  factor  (calculated  based  only  on  grain  orientation 

relative  to  direction  of  applied  uniaxial  stress)  for  each  type  of  slip  system  relates  to  the  extreme  value  volume 
averaged  FS  FIP  for  each  simulation  for  Microstructure  A.  The  other  microstructures  demonstrated  similar  trends. 
In  most  cases,  the  grains  with  the  maximum  grain  averaged  FIP  are  oriented  for  easy  basal  slip  or  in  a  softer 
configuration  relative  to  the  loading  direction.  Additionally,  maximum  Schmid  factors  are  also  observed  for 

prismatic  slip  in  a  few  cases  for  pyramidal  {a^j  and  pyramidal  (a  +  c }  slip.  When  the  magnitude  of  the  Schmid 
factor  for  pyramidal  (a'j  and  pyramidal  (a  +  cj  slip  is  high,  the  magnitude  of  the  Schmid  factor  for  prismatic  slip 
is  also  fairly  high.  This  is  significant  because  prismatic  slip  is  an  easier  slip  mode  than  either  pyramidal  (a'j  and 

pyramidal  (a  +  c^  slip;  therefore,  prismatic  slip  will  likely  be  favored  even  when  the  resolved  shear  stresses  on 
both  slip  systems  is  similar. 

As  noted  previously,  the  primary  a  grain  size  has  been  reported  to  play  a  role  in  the  processes  of  fatigue  crack 
formation.  Figure  13  shows  the  cross  correlation  between  the  grain  volume  associated  with  the  center  of  the 
volume  averaged  extreme  value  FS  FIP  versus  the  extreme  value  FS  FIPs  for  Microstructure  A.  The  volume 
averaged  FIPs  in  this  case  are  calculated  over  the  averaging  volume  of  9.83xl0'5  with  equivalent  grain  size  of 
0.032mm;  the  mean  primary  a  grain  size  for  Microstructure  A  is  0.025mm  with  an  approximate  volume  of 
4.56xl0'5  as  calculated  via  Equation  (5)  (see  Table  1).  The  grains  with  the  extreme  value  FIPs  tend  to  be  larger  than 


15 


Simulated  Microstructure-Sensitive  Extreme  Value  Probabilities  for  High  Cycle  Fatigue  of  Duplex  Ti-6AI-4V 

Craig  P.  Przybyla  and  David  L  McDowell 


the  average  grain  size,  but  the  largest  observed  FIPs  do  not  necessarily  occur  in  the  largest  grains.  This  is  in 
agreement  with  what  has  been  observed  experimentally  in  Ti-6246  [16,  17]. 

The  extreme  value  marked  correlation  functions  are  also  constructed  to  consider  the  crystallographic  attributes 
that  most  strongly  influence  the  predicted  extreme  value  response.  As  demonstrated  previously  [31],  the 
microstructure  dependence  of  the  extreme  value  FIPs  are  characterized  in  terms  of  the  apparent  Schmid  factor  for 
the  different  phases  considered.  Specifically,  Schmid  factors  are  identified  that  correlate  with  the  extreme  value 
response  using  the  previously  introduced  extreme  value  marked  radial  correlation  function  for  the  Schmid  Factor 
(i.e.,  ^max(«)  (msa,mSa  |  r,fl)  )•  By  considering  the  extreme  value  marked  correlation  functions  for  apparent  Schmid 

factors,  the  probabilities  of  specific  grain  orientations  and  misorientations  existing  coincident  with  the  observed 
extreme  value  FIPs  are  indirectly  considered.  To  construct  ^max(«)  ^msa,msa'  |r?Q),  a  tessellated  scale  of  the  Schmid 

factor  in  0.05  increments  was  applied.  As  shown  previously  in  Figure  12,  grains  oriented  for  easy  basal  slip  were 
most  common  at  locations  of  the  estimated  extreme  value  FIPs.  Thus,  to  reduce  the  dataset  to  a  manageable  size, 
the  correlations  between  the  Schmid  factors  were  considered  for  basal  slip  between  0.45  and  0.5  for  the  primary  a 
phase,  and  Schmid  factors  between  0.45  and  0.5  for  the  same  and  for  the  other  types  of  slip  and  phases,  as  shown 
in  Figures  14  and  15  for  Microstructure  A.  The  correlations  that  occurred  with  the  highest  significance  in 
Microstructure  A  are  shown  in  Figure  14a,  which  plots  correlations  between  the  primary  a  grains  oriented  for  easy 
basal  slip  (i.e.,  with  a  Schmid  factor  for  basal  slip  between  0.45  and  0.5)  and  similarly  oriented  grains  of  the  same 
phase.  In  Figure  15d,  the  correlations  are  plotted  that  occurred  with  the  second  highest  probability  in 
Microstructure  A.  These  correlations  describe  the  probabilities  that  primary  a  grains  oriented  for  easy  basal  slip 
occur  near  a+P  colonies  oriented  favorably  for  bcc  slip.  Note  that  the  active  bcc  slip  systems  considered  in  the 
crystal  plasticity  model  are  aligned  with  certain  HCP  slip  systems  according  to  (0001)  //{lOl}  and 


(OOOl)  I/(l  1 1)  in  order  to  maintain  the  BOR.  In  Figures  16-18,  the  correlated  Schmid  factors  are  similarly  plotted 


for  the  different  phases  that  occurred  with  the  highest  probability  in  microstructures  B,  C,  and  D,  respectively. 

In  all  microstructures,  it  was  found  that  primary  a  grains  oriented  with  a  Schmid  factor  for  basal  slip  between 
0.45  and  0.5  (easy  slip)  are  much  more  likely  to  be  found  at  the  location  of  extreme  value  response  than  in  the 
overall  microstructure.  This  agrees  with  what  was  found  by  Bridier  et  al.  [13,  14].  When  considering  multiple 
correlated  microstructure  attributes,  it  was  much  more  probable  to  find  a  primary  a  grain  oriented  for  easy  basal 
slip  near  hard  oriented  a+P  colonies  with  high  Schmid  factors  for  BOR-modified  bcc  slip  in  microstructures  A  and  B 

or  primary  a  grains  with  high  Schmid  factors  for  pyramidal  (a  +  cj  slip  in  microstructures  B  and  D.  The  latter 


previously  mentioned  correlated  orientation  relationships  were  also  identified  by  Gilbert  and  Piehler  [15]  and 
Bantounas  et  al.  [20]  in  Ti-6AI-4V.  Previously  cited  work  [21-26]  demonstrated  that  these  types  of  hard-soft  grain 
combinations  can  dramatically  increase  the  local  stress  state  which  in  turn  promotes  increased  plasticity  in  these 
regions.  It  is  not  surprising  that  simulations  from  microstructures  A  and  C  predict  these  hard-soft  grain  correlations 
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to  occur  with  highest  probability  when  the  primary  a  grains  oriented  for  easy  slip  are  near  the  hard  oriented  a+P 
colonies  when  the  volume  fraction  of  the  primary  a  grains  is  low.  In  microstructures  B  and  D,  the  hard-soft 
correlations  that  occurred  with  the  highest  probably  were  observed  to  be  between  the  primary  a  grains  oriented 
for  slip  near  the  hard  oriented  primary  a  grains,  which  is  expected  when  the  volume  fraction  of  the  primary  a 
grains  is  high  (i.e.,  ~70%).  Other  hard-soft  grain  relationships  in  microstructure  are  also  predicted  to  occur  with 
significantly  higher  probabilities  at  sites  of  extreme  value  response  than  in  the  overall  microstructure.  For  brevity 
only  those  in  Microstructure  A  are  shown  in  Figures  14c)  and  14d),  which  describe  correlations  between  primary  a 

grains  oriented  for  easy  basal  slip  near  other  primary  a  grains  oriented  for  harder  pyramidal  ^ a }  or  pyramidal 

(a  +  cj  slip  modes,  respectively.  Figure  15  shows  correlations  between  primary  a  grains  oriented  for  easy  basal 

slip  and  a+P  colonies  that  in  all  cases  exhibit  harder  slip  modes.  In  short,  the  same  key  correlated  microstructure 
attributes  (in  this  case  grains/phases  with  a  specific  Schmid  factor)  were  identified  across  the  range  of  all  the 
microstructures.  The  only  significant  differences  among  the  different  microstructures  between  the  key  correlated 
attributes  identified  by  the  extreme  value  marked  correlations  functions  were  the  magnitudes  of  the  probabilities 
that  these  correlated  attributes  exist  in  the  regions  marked  by  the  extreme  value  behavior  compared  to  the  overall 
microstructure. 

The  present  study  illustrates  how  this  newly  devised  methodology  can  be  applied  to  a  complex  alloy  system. 
Certain  limitations  are  recognized  in  terms  of  applicability  to  specific  a+p  Ti  alloys.  In  all  cases  the  generated  SVEs 
in  this  study  were  fit  to  specific  grain  size,  orientation  and  disorientation  distributions  that  could  be  very  different 
from  certain  material  systems  examined  experimentally.  For  example,  the  grain  orientation  and  grain 
disorientation  distributions  were  assumed  to  be  random;  however,  this  is  not  very  likely  in  a  real  material  system. 
A  random  ODF  and  DDF  was  chosen  in  this  study  to  facilitate  sampling  of  as  many  different  arrangements  of  grains 
as  possible  to  better  identify  the  most  important  correlated  attributes  related  to  the  processes  of  fatigue  crack 
formation.  In  textured  a+P  Ti  alloys,  some  of  the  arrangements  identified  here  exist  with  little  or  no  probability. 
There  is  no  conceptual  difficulty  in  applying  this  framework  to  compare  effects  of  texture  intensity  and  multiple 
components  of  texture  on  fatigue  response.  Additionally,  in  the  present  study  only  subsurface  material  volumes 
were  considered  since  3D  periodic  boundary  conditions  were  used  on  SVEs.  Surface  effects  on  the  driving  forces 
for  fatigue  crack  formation  have  been  neglected  in  this  work.  In  most  cases,  there  are  competing  surface  and 
subsurface  modes  of  fatigue  crack  formation. 

4.  Conclusions 

This  work  has  demonstrated  the  use  of  a  new  microstructure-sensitive  extreme  value  probabilistic  framework 
that  links  the  extreme  value  driving  forces  for  fatigue  crack  formation  from  polycrystal  plasticity  simulations 
conducted  on  duplex  Ti-6AI-4V  to  the  important  correlated  microstructure  attributes  that  most  influence  that 
extreme  value  response.  Specifically,  it  is  concluded  that: 
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•  The  ellipsoidal  based  microstructure  generator  can  be  a  useful  tool  to  reproduce  synthetic  microstructure 
volumes  that  match  predefined  distributions  of  grain  size,  grain  orientation  and  grain  disorientation  for 
two-phase  polycrystalline  material  systems. 

•  Multiple  statistical  volume  elements  (SVEs)  can  be  constructed  and  simulated  using  appropriate 
constitutive  models  to  estimate  extreme  value  distributions  of  key  response  parameters. 

•  Cyclic  plastic  strain  based  fatigue  indicator  parameters  calculated  over  nonlocal  averaging  volumes  on  the 
scale  of  the  microstructure  attributes  relevant  to  the  mechanisms  of  fatigue  crack  formation  can  be  used 
to  identify  key  microstructure  attributes  that  correlate  with  extreme  value  response. 

•  In  agreement  with  experiments,  primary  a  grains  oriented  favorably  for  basal  slip  (with  a  Schmid  factor 
for  basal  slip  between  0.45  and  0.5)  have  a  much  higher  probability  of  being  associated  with  the  regions 
of  extreme  value  FIPs  than  in  the  overall  microstructure. 

•  In  agreement  with  experiments,  hard-soft  grain  interactions  between  primary  a  grains  oriented  favorably 
for  basal  slip  and  primary  a  grains  oriented  favorably  for  pyramidal  (a  +  c }  or  a+(3  colonies  oriented 

favorably  for  BOR  modified  bcc  slip  correlate  with  increased  driving  forces  for  fatigue  crack  formation 
over  a  relatively  wide  range  of  a+P  Ti  alloy  systems. 
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Tables 

Table  1  Microstructures  for  Ti-6AI-4V  along  with  the  assumed  log-normal  fits  for  the  grain  volume  distributions. 

Assumed  Mean  and  St.  Dev.  for 
Boeing  Supplied  Microstructures  Grain  Size  Distributions 

a+|3  Colony  Primary  a 


Transformed 

Primary  a 

Vol  % 

<7 

a 

Micro.  Name 

|3  Size  (urn) 

Size  (urn) 

Primary  a 

(kirn) 

(kirn) 

(kirn) 

(urn) 

A  Fine  bi-modal  low  a 

50 

10-50 

30% 

50 

5 

25 

10 

B  Fine  bi-modal  high  a 

50 

10-50 

70% 

50 

5 

25 

10 

C  Coarse  bi-modal  low  a 

80 

40-60 

30% 

80 

10 

50 

5 

D  Coarse  bi-modal  high  a 

80 

40-60 

70% 

80 

10 

50 

5 

Table  2  Range  of  number  of  generated  grains  in  cubiodal  SVEs  dimensioned  0.400  mm  along  each  edge  for  the  four 
microstructures. 


# 

Name 

Number  of  Grains 

A 

Fine  bi-modal  low  a 

325-375 

B 

Fine  bi-modal  high  a 

625-675 

C 

Coarse  bi-modal  low  a 

75-100 

D 

Coarse  bi-modal  high  a 

125-175 

Table  3  Size  of  three  different  grain  scale  averaging  volumes  used  to  calculate  FIPs.  The  equivalent  grain  size  is 
calculated  using  Equation  (5). 


Ave.  Vol. 

#  of 

Elements 

Total  Volume 

(mm3) 

Equiv.  Grain  Size 

(mm) 

1 

27 

9.83xl0‘5 

0.032 

2 

125 

4.55xl0‘4 

0.054 

3 

343 

1.25xl0"3 

0.076 
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Table  4  Fitting  parameters  for  Microstructures  A-D  for  the  Gumbel  distribution  for  the  different  averaging  volumes 
defined  in  Table  3. 


Ave.  Vol. 

Micro. 

an 

Un 

K)"Vm» 

R2 

A 

1.257xl07 

9.387x10  8 

0.85 

0.98 

B 

6.562x10s 

1.117x10  7 

1.36 

0.79 

1 

C 

5.032xl06 

1.008xl0"7 

1.97 

0.93 

D 

4.968x10s 

1.485xl0"7 

1.36 

0.88 

A 

4.631xl07 

2.369x10  8 

0.91 

0.97 

B 

1.837xl07 

2.429x10  8 

2.24 

0.67 

2 

C 

1.474xl07 

3.304x10  8 

2.05 

0.92 

D 

1.282xl07 

4.401x10  8 

1.77 

0.82 

A 

1.221xl08 

9.443x10  9 

0.87 

0.97 

B 

4.848xl07 

1.046x10  8 

1.97 

0.68 

3 

C 

3.734x10 

1.342x10 

2.00 

0.92 

D 

3.004xl07 

1.760x10  8 

1.89 

0.82 

Table  5  Fitting  parameters  for  Microstructures  A-D  for  the  Frechet  distribution. 


Ave.  Vol. 

Micro. 

k 

R2 

A 

1.634E+00 

7.719E-08 

0.91 

1 

B 

1.998E+00 

1.212E-07 

0.97 

C 

9.570E+00 

1.907E-01 

0.89 

D 

1.599E+00 

1.378E-07 

0.96 
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Figures 


Cyclic  Plasticity  and  Ratcheting 


> 

V) 

< Si 

S 

H 

55 

L 

— ► 

i 

H 

E 

Elastic  Cyclic  Ratcheting 

Shakedown  Plasticity 


Figure  1:  Schematic  showing  definitions  of  reversed  cyclic  plasticity  ( Asp  )  ^  and  ratcheting  ( A Sp  ^ 


/ratch 
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a) 


b) 


Figure  2:  a)  A  2-D  schematic  demonstrating  the  packing  of  ellipses  and  b)  a  2-D  section  of  a  microstructure 
generated  using  the  ellipsoid  packing  algorithm  with  a  bi-modal  grain  size  distribution. 
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Figure  3:  Aspect  ratios  of  the  initial  generated  ellipsoids  (dots)  and  grains  (circles)  after  uniform  grain  growth. 
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Figure  4:  Target  log-normal  and  fit  model  distributions  of  the  grain  size  for  the  a+(3  colonies  (pj  and  primary  a 
grains  (p2)  for  a  randomly  generated  fine  bi-modal  microstructure. 
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Figure  5:  Target  log-normal  and  fit  model  distributions  of  the  grain  size  for  the  a+(3  colonies  (pi)  and  primary  a 
grains  (p2)  for  a  randomly  generated  coarse  bi-modal  microstructure. 
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Figure  6:  Meshed  cuboidal  SVE  with  an  edge  dimension  of  0.400  mm  for  an  arbitrary  instantiation  of 
Microstructure  A  with  contours  of  the  Mises  stress  at  the  peak  strain  of  0.6%  after  10  cycles. 
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10 

N  Cycles 


15 


20 


c) 


Figure  7:  The  percentage  of  primary  a  grain  and  a+(3  colonies  that  exhibit  an  effective  plastic  strain  above  a) 
l.OxlO"10,  b)  l.OxlO"8,  and  c)  l.OxlO"6  over  for  the  first  20  cycles  at  0.6%  maximum  strain  for  R  =  0  cycling. 
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FSFIP 


c) 


Figure  8:  The  extreme  value  FIP  calculated  over  a  cube  shaped  averaging  volume  with  equivalent  grain  size  of 
0.032mm  for  Microstructures  A-D  plotted  on  a  Gumbel  probability  scale  for  a)  the  cumulative  effective  plastic 
strain  (CEPS),  b)  the  maximum  plastic  shear  strain  (MPSS),  and  c)  the  Fatemi-Socie  (FS)  FIPs. 


33 


Simulated  Microstructure-Sensitive  Extreme  Value  Probabilities  for  High  Cycle  Fatigue  of  Duplex  Ti-6AI-4V 

Craig  P.  Przybyla  and  David  L.  McDowell 


FSFIP 


a) 


b) 


Figure  9:  The  extreme  value  Fatemi-Socie  (FS)  FIP  calculated  over  a  cube  shaped  averaging  volume  with  equivalent 
grain  size  of  0.032mm  for  the  four  analyzed  microstructures  plotted  on  a)  Gumbel  (Type  I),  b)  Frechet  (Type  II) 
probability  scales.  Note  that  the  equations  given  for  the  least  squares  linear  regression  are  such  that 
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y  =  ln(l/ln(l/ p)^)  and  x  =  FIP  for  the  Gumbel  probability  plot  and  x  =  ln(FIP)  for  the  Frechet  probability  plot, 
respectively. 
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Figure  10  The  extreme  value  Fatemi-Socie  (FS)  FIP  calculated  over  a  cube  shaped  averaging  volume  with  equivalent 
grain  size  of  0.032mm  for  Microstructure  A  on  a  Gumbel  probability  scale  for  the  first  10  cycles.  Note  that  the 

equations  given  for  the  least  squares  linear  regression  are  such  that  y  =  ln(l/ln(l/ andx  =  FIP. 
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FSFIP 


a) 


FSFIP 


b) 


Figure  11:  The  extreme  value  Fatemi-Socie  (FS)  FIP  calculated  over  averaging  volume  with  equivalent  grain  size  of 
a)  0.054mm  and  b)  0.076mm  for  the  range  of  microstructure  considered  on  a  Gumbel  probability  scale.  Note  that 
the  equations  given  for  the  least  squares  linear  regression  are  such  that  y  =  ln(l/ln(l/  and  x  =  FIP. 
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Figure  12:  The  Schmid  factors  (calculated  based  on  the  geometry  of  slip  system  relative  to  the  loading  direction) 
for  basal,  prismatic,  pyramidal  <a>,  and  pyramidal  <a+c>  for  each  extreme  value  Fatemi-Socie  FIP  calculated  over  a 
cube  shaped  averaging  volume  with  equivalent  grain  size  of  0.032mm  estimated  from  the  100  simulated  SVEs  for 
Microstructure  A. 
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Figure  13:  The  extreme  value  Fatemi-Socie  FIPs  calculated  over  a  cube  shaped  averaging  volume  with  equivalent 
grain  size  of  0.032mm  plotted  against  the  grain  volume  of  the  grain  coincident  with  the  extreme  value  FIP  for  the 
100  simulated  SVEs  for  Microstructure  A.  The  average  grain  volume  for  Microstructure  A  calculated  via  Equation 

(5)  based  on  the  average  grain  size  of  0.025mm  is  shown  by  V  . 
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Figure  14:  Radial  correlation  functions  for  the  complete  microstructure  compared  with  extreme  value  marked 
radial  correlation  functions  describing  the  correlation  between  Schmid  factors  ma  for  basal  slip  between  0.45  and 

0.5  for  the  primary  a  phase  and  the  Schmid  factors  rria  for:  a)  basal  slip  between  0.45  and  0.5,  b)  prismatic  slip 
between  0.45  and  0.5,  c)  pyramidal  ^ a ^  slip  between  0.45  and  0.5,  and  d)  pyramidal  (a  +  c ^  slip  between  0.45 
and  0.5  for  the  primary  a  phase  in  Microstructure  A. 
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Figure  15  Radial  correlation  functions  for  the  complete  microstructure  compared  with  extreme  value  marked  radial 
correlation  functions  describing  the  correlation  between  the  apparent  Schmid  factors  ma  for  basal  slip  between 

0.45  and  0.5  for  the  primary  a  phase  and  the  apparent  Schmid  factors  ma  for:  a)  basal  slip  between  0.45  and  0.5, 

b)  prismatic  slip  between  0.45  and  0.5,  c)  pyramidal  ^ a }  slip  between  0.45  and  0.5,  and  d)  (\  1  l^jl  10}  bcc  slip 

(transformed  into  the  hexagonal  coordinate  system  via  the  BOR)  between  0.45  and  0.5  for  the  a+(3  colony  phase  in 
Microstructure  A. 
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Figure  16  Radial  correlation  functions  for  the  complete  microstructure  compared  with  extreme  value  marked  radial 
correlation  functions  describing  the  correlation  between  the  apparent  Schmid  factors  ma  for  basal  slip  between 

0.45  and  0.5  for  the  primary  a  phase  and  the  apparent  Schmid  factors  ma  for:  a)  basal  slip  between  0.45  and  0.5 
and  b)  pyramidal  (a  +  c }  slip  between  0.45  and  0.5  for  the  primary  a  phase  in  Microstructure  B. 
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a) 


b) 


Figure  17  Radial  correlation  functions  for  the  complete  microstructure  compared  with  extreme  value  marked  radial 
correlation  functions  describing  the  correlation  between  the  apparent  Schmid  factors  ma  for  basal  slip  between 

0.45  and  0.5  for  the  primary  a  phase  and  the  apparent  Schmid  factors  ma  for:  a)  basal  slip  between  0.45  and  0.5 
and  b)  l  l^{l  10}  bcc  slip  (transformed  into  the  hexagonal  coordinate  system  via  the  BOR)  between  0.45  and  0.5 
for  the  a+P  colony  phase  in  Microstructure  C. 
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a) 


b) 


Figure  18  Radial  correlation  functions  for  the  complete  microstructure  compared  with  extreme  value  marked 
radial  correlation  functions  describing  the  correlation  between  the  apparent  Schmid  factors  ma  for  basal  slip 

between  0.45  and  0.5  for  the  primary  a  phase  and  the  apparent  Schmid  factors  ma  for:  a)  basal  slip  between  0.45 
and  0.5  and  c)  pyramidal  (a  +  c ^  slip  between  0.45  and  0.5  for  the  primary  a  phase  in  Microstructure  D. 
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